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Abstract 


In  this  report  we  describe  experiments  carried  out  at  the  DRDC  Atlantic 
tank  facility  with  the  high  frequency  broadband  Multi-Mode  Pipe  Projector 
(MMPP).  In  particular,  we  examine  the  fidelity  with  which  one  can  produce 
specified  broadband  pulses  with  and  without  compensation  for  the  transducer’s 
spectral  characteristics.  It  is  shown  that  extended  pulses  or  pulse  sequences 
can  be  accurately  produced  and  that  matched  filtering  very  effectively  com¬ 
presses  the  time  series  in  the  case  of  extended  pulses.  In  addition,  the  vertical 
beampattern  of  the  transducers  is  discussed  and  the  results  of  some  experi¬ 
ments  with  2  transducers  are  also  presented. 


Resume 


Le  present  rapport  comprend  la  description  d’experiences  effectuees  dans  le 
reservoir  acoustique  de  RDDC  Atlantique  en  utilisant  le  projecteur  a  tube 
multimode  (MMPP  pour  Multi-Mode  Pipe  Projector)  haute  frequence  a  large 
bande.  En  particulier,  la  fidelite  avec  laquelle  on  peut  reproduire  des  im¬ 
pulsions  a  large  bande  specifiees  en  compensant,  et  sans  compenser,  pour  les 
caracteristiques  spectrales  du  transducteur  est  examinee.  Les  impulsions  sont 
generees,  et  le  filtrage  adapte  permet  de  comprimer  de  facon  tres  efficace  la 
serie  chronologique  des  impulsions  de  duree  prolongee.  Le  rapport  comporte 
egalement  une  discussion  du  diagramme  de  faisceau  vertical  des  transducteurs 
et  les  resultats  de  certaines  experiences  elfect  uees  avec  deux  transducteurs. 
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Executive  summary 


BACKGROUND 

This  is  the  first  report  in  a  series  of  reports  dealing  with  aspects  of  the  TIF 
project  “Exploiting  Ultra-Wideband  and  Coded  Sonar  Pulses”.  In  this  report 
we  investigate  the  characteristics  and  potential  of  the  smallest  of  the  Multi- 
Mode  Pipe  Projector  (MMPP)  transducers  which  has  a  useful  frequency  range 
from  about  20-120  kHz.  This  very  wide  bandwidth  allows  the  transducer  to 
accurately  create  either  very  short  duration  pulses  or  complicated  extended 
pulses  such  as  Chirps  or  sequences  of  smaller  pulses.  It  is  hoped  that  the 
wideband  capabilities  of  the  MMPP  technology  will  be  useful  in  a  variety  of 
applications  such  as  underwater  communications,  decoy  signals,  low  amplitude 
extended  sonar  pulses,  and  target  scattering.  These  applications  have  naval 
relevance  for  improved  mine  classification,  stealthy  underwater  communica¬ 
tion,  torpedo  jamming,  and  diver  detection  systems. 

SIGNIFICANCE  OF  RESULTS 

It  is  shown  that  the  high  frequency  MMPP  transducer  has  a  wide  frequency 
bandwidth.  The  spectral  characteristics  of  the  transducer  can  be  compensated 
for  in  the  construction  of  the  waveform  input  to  the  transducer.  The  resulting 
output  waveform  accurately  replicates  the  desired  waveform.  Matched  filter¬ 
ing  techniques  can  be  effectively  applied  to  the  time  series  recorded  from  the 
transducer  to  increase  the  signal  to  noise  ratio  and  to  temporally  compress  the 
received  waveform.  The  beampattern  (horizontal  and  vertical)  characteristics 
of  the  transducer  vary  significantly  with  the  orientation  of  the  transducer  and 
for  different  applications,  it  may  be  advantageous  to  vary  the  orientation  of  the 
transducer.  The  characteristics  of  the  transducer,  as  described  in  this  report, 
indicate  that  it  may  be  useful  for  a  variety  of  important  naval  applications 
such  as  underwater  communications,  torpedo  decoy  signals,  diver  detection 
systems,  and  mine  classification. 

FUTURE  WORK 

There  is  much  work  remaining  with  the  use  of  the  MMPP  transducers.  The 
characteristics  of  the  other  2  sized  transducers  will  be  investigated.  As  well,  the 
use  of  these  transducers  in  conjunction  with  each  other  to  increase  and  flatten 
the  overall  spectral  response  will  also  be  investigated.  The  application  of 
these  transducers  and  possibly  arrays  of  these  transducers  will  be  investigated 
in  more  realistic,  longer  range  sea  trials. 

Fawcett,  J.,  Sildam,  .1.,  Miller,  T.,  Fleming,  R.,  and  Trevorrow,  M.,  2005. 
Broadband  pulse  synthesis  with  the  high  frequency  Multi-Mode  Pipe  Projec¬ 
tor.  DRDC  Atlantic  TAI  2004-272,  Defence  R&D  Canada  -  Atlantic. 
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Sommaire 


INTRODUCTION 

Le  present  rapport  est  le  premier  d’une  serie  de  rapports  qui  vont  traiter  de 
divers  aspects  du  projet  FIT  portant  sur  l’utilisation  d’impulsions  sonar  a 
bande  ultra-large  et  des  impulsions  sonar  codees  (Exploiting  Ultra- Wideband 
and  Coded  Sonar  Pulses).  Des  travaux  de  recherche  sur  les  caracteristiques  et 
le  potentiel  du  plus  petit  des  transducteurs  MMPP,  dont  la  plage  de  frequences 
utiles  s’etend  d’environ  20  a  120  kHz,  sont  decrits  dans  le  present  rapport. 
Cette  largeur  de  bande  etendue  permet  au  transducteur  de  creer  avec  precision 
soit  des  impulsions  de  tres  courte  duree,  soit  des  impulsions  complexes  de  duree 
prolongee,  conime  des  ondes  entretenues  modulees  en  frequence  (chirps)  ou  des 
sequences  d’impulsions  plus  courtes.  On  espere  que  les  capacites  large  bande  de 
la  technologie  MMPP  seront  utiles  pour  de  nombreuses  applications,  dont  les 
telecommunications  sous-marines,  les  signaux  de  leurre,  les  impulsions  sonar 
de  duree  prolongee  a  faible  amplitude,  et  la  dispersion  par  la  cible.  Ces  appli¬ 
cations  presentent  un  interet  pour  la  marine  en  ce  qui  a  trait  a  la  classification 
amelioree  des  mines,  aux  telecommunications  sous-marines  indetectables,  aux 
systemes  de  brouillage  visant  les  torpilles,  et  aux  systemes  de  detection  de 
plongeurs. 

PORTEE 

II  est  demontre  que  le  transducteur  MMPP  haute  frequence  a  une  grande 
largeur  de  bande.  On  peut  compenser  pour  les  caracteristiques  spectrales  du 
transducteur  par  la  generation  de  formes  d’onde  d’entree  au  transducteur. 
La  forme  d’onde  de  sortie  resultante  correspond  avec  precision  a  la  forme 
d’onde  desired.  Des  techniques  de  filtrage  adapte  peuvent  etre  appliquees  ef- 
ficacement  a  la  serie  chronologique  enregistree  a  partir  du  transducteur  afin 
d’augmenter  le  rapport  signal  sur  bruit  et  de  comprimer  temporairement  la 
forme  d’onde  recue.  Les  caracteristiques  des  diagraming  de  faisceau  (hori¬ 
zontal  et  vertical)  du  transducteur  varient  de  facon  significative  en  fonction 
de  son  orientation.  Pour  diflerentes  applications,  il  peut  done  etre  avantageux 
de  modifier  l’orientation  du  transducteur.  Les  caracteristiques  du  transduc¬ 
teur,  decrites  dans  le  present  rapport,  indiquent  que  cela  pourrait  etre  utile 
dans  le  cas  de  nombreuses  applications  importantes  pour  la  marine,  comme 
les  telecommunications  sous-marines,  les  signaux  de  leurre  visant  les  torpilles, 
les  systemes  de  detection  de  plongeurs  et  la  classification  des  mines. 

RECHERCHES  FUTURES 

II  reste  encore  beaucoup  de  travail  a  faire  en  ce  qui  a  trait  a  l’utilisation  des 
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transducteurs  MMPP.  On  compte  faire  des  recherches  sur  les  caracteristiques 
des  deux  autres  transducteurs,  ainsi  que  sur  l’utilisation  combinee  de  ces  trans¬ 
ducteurs  afin  d’augmenter  et  d’aplatir  la  reponse  spectrale  globale.  De  plus, 
l’utilisation  de  ces  transducteurs  et,  possiblement,  de  reseaux  de  transducteurs 
fera  l’objet  d’essais  en  mer  sur  des  plus  grandes  distances,  dans  un  environ- 
nement  plus  realiste. 
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1  INTRODUCTION 


In  2003  TIF  funding  was  granted  for  the  project  “Exploiting  Ultra-Wide-Band 
and  Coded  Sonar  Pulses”.  This  project  proposed  using  the  transducer  technol¬ 
ogy,  Multi-Mode  Pipe  Projector  (MMPP),  developed  at  DRDC  Atlantic  [1,2] 
as  the  type  of  transducer  for  the  study  of  various  wideband  applications.  In 
particular,  3  different  MMPP  sizes  will  be  considered  during  this  TIF  project. 
For  the  studies  in  this  report,  we  focus  on  the  smallest  transducer  which  has 
significant  output  in  the  frequency  range  20-120  kHz.  The  dimensions  of  the 
transducer  are  approximately,  5.3  cm  high  and  a  diameter  of  3.7  cm.  The  small 
size  of  this  transducer  means  that  it  should  be  straightforward  to  construct 
systems  with  arrays  of  these  transducers. 

The  concept  of  the  MMPP  is  based  upon  synthesizing  a  wideband  output  for  a 
transducer  by  constructing  a  transducer  with  a  set  of  many  resonances  spread 
across  the  frequency  spectrum.  This  results  in  a  significant  power  output  over 
the  entire  band  of  frequencies  spanned  by  the  resonances  of  the  system.  It 
also  means,  however,  that  the  resulting  power  spectrum  for  the  transducer  is 
not  flat  as  there  are  peaks  and  nulls  resulting  from  the  combined  response 
of  the  resonances  of  the  system.  Thus,  it  is  not  clear  with  what  fidelity  one 
can  reproduce  specific  pulse  types.  For  some  sonar  applications,  it  may  not 
be  very  important  to  be  able  to  reproduce  exactly  specified  pulses  -  it  may 
be  sufficient  to  simply  use  the  pulse  output  by  the  transducer.  However,  for 
applications  such  as  producing  realistic  decoy  signals,  it  is  clearly  important 
to  be  able  to  accurately  reproduce  specified  signals. 

In  the  first  part  of  Section  2.2  we  show  the  transducer  response  for  constant 
bandwidth  signals  (Sine  functions  in  the  time  domain)  as  a  function  of  the 
centre  frequency  and  bandwidth  for  a  single  transducer.  We  also  show  the 
responses  for  2  transducers  mounted  vertically  and  wired  in  series,  with  re¬ 
versed  and  the  same  polarities  and  wired  in  parallel  with  reversed  and  the 
same  polarity. 

With  a  knowledge  of  the  power  spectrum  output  of  the  transducer  over  a  given 
bandwidth,  we  can  attempt  to  design  a  compensated  input  pulse  which  yields 
the  desired  output  pulse.  The  results  of  this  process  will  be  presented  for 
both  narrow  and  extended  pulses.  As  well,  the  result  of  match-filtering  the 
recorded  time  series  will  be  presented.  Finally,  the  vertical  beampattern  of 
the  transducer  will  be  investigated  and  the  significance  of  the  beampattern 
discussed. 
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2  PULSE  EXPERIMENTS 


2.1  Pulse  types 


In  this  section  we  describe  the  main  pulse  types  we  use  in  the  tank  experiments. 
First,  we  consider  a  constant  bandwidth  pulse  with  centre  frequency  fc  and 
bandwidth  Bf.  In  the  time  domain  the  pulse  is  given  by 


S8(t) 


sin(2i:Bft/2) 

nBft 


cos(27r/cf) 


(1) 


For  very  small  bandwidths  the  signal  is  harmonic  in  nature  with  frequency 
fc,  for  larger  bandwidths  the  pulses  becomes  shorter  and  contain  only  a  few 
cycles  of  the  harmonic  signal.  A  similar  signal  is  obtained  by  considering  a 
Gaussian  distribution  in  the  frequency  domain  (and  its  image  for  the  negative 
frequencies) , 

Sg(t)  =  cos(27r/cf)  exp(— 2ir2ajt2)  (2) 

where  oy  is  the  standard  deviation  of  the  pulse  in  the  frequency  domain.  For 
extended  pulses  we  use  the  CHIRP  subroutine  of  the  signal  processing  toolbox 
in  MATLAB.  Here  the  instantaneous  frequency  of  the  pulse  will  vary  linearly 
from  a  specified  /0  at  t  =  0  to  fi  at  a  specified  end  time  (we  will  use  one 
millisecond) . 


2.2  Frequency  response  measurements 

The  experiments  we  describe  in  this  report  took  place  in  the  DRDC  Atlantic 
acoustic  calibration  tank.  The  high-frequency  MMPP  was  deployed  at  a  depth 
of  1.83  m  and  a  B&K  8100  calibrated  hydrophone  was  deployed  2  m  away.  In 
Fig.l  we  show  the  Transmit  Voltage  Response  curves  obtained  for  the  5  dif¬ 
ferent  high  frequency  MMPP  transducers  using  the  standard  DRDC  Atlantic 
calibration  acquisition  and  software  suite.  As  can  be  seen  the  5  curves  are  very 
similar  except  for  a  deeper  null  in  one  curve  in  the  region  of  80  kHz.  There  are 
significant  vertical  beampattern  effects  at  the  high  frequencies  and  the  differ¬ 
ence  here  may  be  due  to  a  slight  change  in  the  geometry  or  relative  orientation 
of  the  transducer.  However,  there  is  certainly  a  null  in  the  responses  of  all  5 
transducers  at  approximately  30  and  80  kHz. 

A  data  acquisition  and  waveform  generation  program  was  written  in  LAB- 
VIEW  which  allowed  for  the  sequential  generation  of  Sine  and  Gaussian  pulses 
as  a  function  of  centre  frequency  and  bandwidth  (with  a  specified  oj  in  the 
case  of  the  Gaussian  pulses)  and  for  the  recording  of  the  resulting  signals  at 
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Figure  1:  Transmit  voltage  response  (dB  re  1  p  Pa/Volt  @  lm)  curves  for  the 
MMPP  transducer  at  1.83  m  (approx.  6  feet)  depth 


the  hydrophone.  For  the  data  of  this  report,  the  time  series  was  sampled  at 
250000  samples/second.  As  well,  text  files  of  arbitrary  waveforms  can  be  used 
by  the  LAB  VIEW  program. 

In  order  to  test  transducer  linearity,  a  series  of  tests  with  wideband  pulses 
were  conducted.  In  Fig.  2  we  show  a  two-dimensional  image  of  the  resulting 
power  spectra,  computed  by  a  MATLAB  program,  as  the  centre  frequency  fc 
of  a  Sine  pulse  is  varied  from  10  to  100  kHz  in  steps  of  200  Hz  for  a  bandwidth 
of  80  kHz.  The  x-coordinate  is  the  centre  frequency  of  the  Sine  pulse  in 
kHz  and  the  z-coordinate  is  the  frequency  axis  for  the  computed  FFT-power 
spectra.  As  can  be  seen,  the  input  Sine  functions  do  not  excite  energy  in  the 
band  outside  the  80  kHz  band  width.  Another  encouraging  feature  is  that 
at  a  given  frequency  value  for  the  spectral  response,  the  value  is  essentially 
independent  of  the  centre  frequency  fc  (as  long  as  it  lies  within  the  ±  40  kHz 
of  the  centre  frequency).  A  somewhat  unusual  feature  is  the  apparent  folding 
of  energy  at  the  upper  left  and  lower  right  corners  of  the  plot.  For  the  case  of 
the  low  frequencies  this  is  due  to  the  fact  that  some  of  the  frequency  response 
centred  about  the  negative  value  of  the  centre  frequency  appears  in  the  positive 
spectrum.  The  effect  in  the  lower  right  is  due  to  spectral  aliasing.  In  Fig.  3  the 
spectrum  for  a  centre  frequency  of  50  kHz  and  total  bandwidth  of  80  kHz  is 
shown.  As  can  be  seen  this  spectrum  is  very  similar  to  the  calibration  spectra 
of  Fig.  1.  Thus,  it  seems  that  the  transducer  is  exhibiting  predictable  linear 
behaviour. 
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Figure  2:  The  spectral  response  of  the  MMPP  (HF)  transducer  (z-axis)  as  a 
function  of  the  centre  frequency  of  the  80-kHz  constant  bandwidth  signal 


Due  to  the  small  size  of  the  high-frequency  MMPP,  it  is  straightforward  to 
form  arrays  of  these  transducers.  In  Fig.  4  a  2  element  vertical  array  with  a 
rigid  connection  of  approximately  4  cm  spacing  is  shown.  These  2  transducers 
can  be  wired  in  series  or  in  parallel  and  they  can  transmit  in  or  out  of  phase 
with  each  other.  In  Fig.  5  we  show  the  two-dimensional  spectral  response  for 
the  4  possible  arrangements.  In  this  plot,  the  top  row  are  the  two  transducers 
wired  in  series,  out-of-phase  and  in-phase.  The  second  row  shows  the  results 
for  the  transducers  wired  in  parallel.  It  is  important  to  note  that  these  results 
have  all  been  normalized  by  their  maximum  amplitude  and  thus  these  plots 
do  not  indicate  the  absolute  levels  of  the  system  output.  In  Figs. 6  and  7  we 
compare  the  spectra  from  the  single  transducer  (blue)  and  the  two-transducer 
system  (parallel,  out-of-phase)  (red)  for  centre  frequencies  of  40  and  80  kHz. 

As  can  be  seen,  the  use  of  2  transducers,  particularly  when  used  out-of-phase 
has  a  noticeable  effect.  In  Figs. 6  and  7  the  resulting  frequency  response  for  the 
two-transducer  system  is  flatter  at  the  high  frequencies  of  100  kHz  and  beyond 
and  also  in  Fig. 7  it  can  be  seen  that  the  null  at  30  kHz  has  been  somewhat 
shifted  to  a  lower  frequency.  Thus  it  maybe  possible  to  shape  the  spectral 
response  by  using  more  than  one  of  these  transducers  together.  It  might  be 
useful  in  the  future  to  investigate  more  thoroughly  the  effects  of  combining  the 
individual  transducers  and  the  effect  of  the  spacing  distance  between  them. 
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Figure  3:  The  spectral  response  of  the  MMPP  (HF)  transducer  for  the  centre 
frequency  of  50  kHz  and  bandwidth  of  80  kHz. 


Figure  4:  Two  MMPP  transducers  vertically  aligned 
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Figure  5:  The  spectral  responses  of  2  MMPP  (HF)  transducers  as  a  function  of 
the  centre  frequency  of  the  80-kHz  constant  bandwidth  signal  for  [top  row] 
transducers  wired  in  series,  out-of-phase,  in-phase  [bottom  row]  in  parallel, 
out-of-phase,  in-phase 


Frequency  (kHz) 

Figure  6:  The  spectral  responses  of  a  single  MMPP  (HF)  transducer  (blue)  and 
2-transducer  system  (parallel,  out-of-phase)  (red)  for  the  centre  frequency  of  40 
kHz  and  bandwidth  of  80  kHz. 
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Figure  7:  The  spectral  responses  of  a  single  MMPP  (HF)  transducer  (blue)  and 
2-transducer  system  (parallel,  out-of-phase)  (red)  for  the  centre  frequency  of  80 
kHz  and  bandwidth  of  80  kHz. 

2.3  Uncompensated  and  compensated  time 
series 

In  this  section  we  first  show  in  Fig. 8  some  of  the  spectra  and  time  series  for 
20  kHz  bandwidth  Sine  pulses  with  centre  frequencies  of  30,  50,  70,  and  90 
kHz.  As  well  we  superimpose  the  desired  Sine  pulse  on  the  experimental  time 
series.  The  location  of  the  theoretical  pulse  is  determined  by  finding  the  time 
shift  which  maximizes  the  cross-correlation  of  the  theoretical  pulses  with  the 
recorded  signal.  There  is  no  physical  interpretation  as  to  this  location;  it  is 
simply  the  location  which  mathematically  produces  the  largest  correlation.  As 
is  seen  in  Fig.  8  the  resulting  matches  between  the  theoretical  and  the  recorded 
time  series  are  not,  in  general,  good  for  these  uncompensated  input  waveforms. 
It  is  clear  from  the  results  of  Fig.  8  that  the  output  spectra  resulting  from  a 
constant  input  spectra  (across  the  specified  bandwidth)  are  far  from  flat.  This 
is  not  surprising  considering  the  spectral  response  of  the  transducer  as  shown 
in  Fig.  1.  In  the  top  spectral  plot  of  Fig.  8  there  is  no  apparent  energy  at 
30  kHz  even  though  the  Sine  pulse  is  centred  at  30  kHz.  This  is  due  to  the 
fact,  as  seen  in  Fig.  7,  that  there  is  a  significant  null  at  approximately  30  kHz. 
From  the  relative  flatness  of  the  spectrum,  we  would  expect  to  obtain  the  best 
uncompensated  pulses  for  centre  frequencies  in  the  range  of  approximately  40- 
GO  kHz  and  indeed  the  response  at  50  kHz,  the  second  plot  of  Fig.  8,  is  the 
closest  to  the  desired  output. 


Frequency  (kHz) 
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Figure  8:  The  resulting  power  spectrum  (dB)  for  20-kHz  bandwidth  pulses 
centred,  top  to  bottom,  at  30,50,70,  and  90  kHz  respectively  and  the 
corresponding  time  series  -  red  is  the  original  input  signal.  Both  the  time  series 
and  spectra  have  been  normalized  to  have  a  peak  amplitude  of  unity 

By  considering  the  complex- valued  (i.e. ,  both  amplitude  and  phase)  spectral 
response  of  the  transducer  to  a  given  input  waveform,  we  can  compute  the 
transfer  function  to  apply  to  a  waveform  in  order  to  produce  that  waveform  at 
the  output  of  the  transducer.  Letting  S(f)  be  the  recorded  output  spectrum, 
T(f)  the  transfer  function  of  the  transducer,  and  B(f)  the  original  input 
spectrum,  then  one  can  write 


S{f)  =  T(f)B(f) 


(3) 


Let  us  define 


qu)  =  i  mi)  =  §jjj.  (4) 

From  Eq.(3),  it  is  easy  to  show  that  if  the  waveform  corresponding  to  the 
spectrum  Q(f)B(f)  is  input  then  the  waveform  output  from  the  transducer 
will  correspond  to  B(f).  In  practice,  as  S(f)  becomes  relatively  small  in 
amplitude  (for  example,  with  respect  to  the  background  noise  spectrum)  the 
division  in  Eq.(4)  may  become  unstable.  Thus  we  use  the  following  formula 
for  Q(f)  after  S(f)  has  been  normalized  to  have  a  maximum  amplitude  of 
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Figure  9:  The  recorded  time  series  for  an  uncompensated  Sine  pulse  with  a 
centre  frequency  of  75  kHz  and  a  bandwidth  of  80  kHz. 


unity, 


Q(f) 


B(.f) 

S(/)/|S(/)|(|S(/)|  +  a) 


(5) 


where  a  is  a  scalar  which  we  set  to  0.05  in  the  examples  below.  The  com¬ 
pensating  spectrum  Q(f)  can  be  used  for  any  desired  spectrum  B(f)  which  is 
contained  within  the  same  frequency  bandwidth  as  B(f).  It  should  be  noted 
that  in  this  report  we  use  a  spectral  compensation  approach.  It  is  possible  to 
construct  a  finite  time  series  from  Q(f)  and  perform  the  compensation  in  the 
time  domain  by  a  convolution.  For  some  applications  this  could  be  preferable. 


As  an  example  of  the  spectral  compensation  of  an  input  waveform,  let  us 
consider  the  80-kHz  bandwidth  response  to  the  pulse  with  a  centre  frequency 
of  75  kHz.  The  recorded  resulting  pulse  is  shown  in  Fig.  9  and  its  spectrum  in 
Fig.  10.  It  is  clear  from  Fig. 9  that  the  output  waveform  is  far  from  the  very 
sharp  pulse  which  should  result  from  this  bandwidth  and  it  is  clear  from  Fig.  10 
that  the  spectral  response  is  far  from  fiat.  Performing  the  operation  of  Eq.(5) 
we  obtain  the  results  shown  in  Fig.  11  for  the  amplitude  of  the  compensation 
function  Q(f)  and  in  Fig.  12  for  the  phase.  For  the  phase  of  the  spectrum 
of  Q{f),  we  first  unwrapped  the  phase  of  the  spectrum  (using  the  MATLAB 
function  UNWRAP)  and  then  differenced  the  phase  with  respect  to  frequency 
to  make  the  linear  phase  change  (which  can  result  simply  from  propagation  or 
a  time  delay)  a  constant.  As  can  be  seen,  the  significant  phase  changes  occur 
at  approximately  60  and  80-kHz. 
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Figure  10:  The  experimental  spectrum  for  an  uncompensated  Sine  pulse  with  a 
centre  frequency  of  75  kHz  and  a  bandwidth  of  80  kHz. 


It  is  not  clear  what  the  resulting  power  loss  to  the  output  signal  will  be  as  a 
result  of  the  spectral  flattening.  In  the  normalization  of  Eq.(5)  the  peak  value 
of  S(f)  is  normalized  to  have  unit  amplitude  so  that  the  division  effectively 
emphasizes  the  spectral  components  for  frequencies  other  than  the  peak  fre¬ 
quency.  However,  the  time  series  which  is  supplied  to  the  MMPP  waveform 
generator  is  first  normalized  to  have  a  peak  amplitude  of  unity  so  that  the 
absolute  levels  of  the  spectrum  do  not  affect  the  compensated  time  series  that 
is  computed.  However,  by  decreasing  the  relative  contribution  of  frequencies 
(e.g.,  40  kHz)  and  increasing  the  spectral  amplitude  of  “inefficient”  frequen¬ 
cies  (e.g.,  80  kHz)  one  would  expect  the  transducer  to  produce  less  acoustic 
energy  for  a  fixed  amount  of  applied  voltage.  For  example,  the  difference  be¬ 
tween  the  response  at  42  kHz  to  that  at  approximately  78  kHz  is  about  15 
dB.  The  linear  amplitude  reduction  for  the  42  kHz  response  corresponding  to 
a  7.5  dB  drop  (i.e. ,  the  output  at  42  kHz  is  decreased  by  7.5  dB  and  the  ouput 
at  78  kHz  increased  by  the  same  amount)  is  0.421.  The  overall  output  levels 
can  be  increased  by  increasing  the  voltage  applied  to  the  system  and  thus,  in 
practice,  the  levels  which  can  be  obtained  with  the  spectrally-flattened  signal 
will  depend  upon  the  maximum  voltage  levels  which  can  be  applied  to  the 
transducer. 

In  Fig.  13  we  show  the  compensated  signal  used  by  the  MMPP  for  the  80-kHz 
Sine  pulse  (35  kHz  <  /  <  115  kHz).  It  is  clearly  significantly  different  than 
the  desired  Sine  pulse.  The  resulting  recorded  signal  (and  the  overlaid  desired 
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Figure  11:  The  spectral  amplitude  of  the  compensating  spectrum  Q(f) 


Figure  12:  The  phase  of  the  compensating  spectrum  Q(f) 
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Figure  13:  The  pre-compensated  time  series  input  to  the  waveform  generator  to 
produce  a  80-kHz  Sine  Function 


Sine  pulse)  is  shown  in  Fig.  14  and  its  spectrum  in  Fig.  15.  From  Fig.  14  it  can 
be  seen  that  the  compensation  was  largely  successful  in  producing  the  desired 
output  waveform.  Fig.  15  shows  that  the  resulting  spectrum  is  quite  flat  with 
the  exception  of  frequencies  greater  than  approximately  100  kHz  where  the 
relative  levels  are  somewhat  low.  The  lower  levels  are  probably  due  to  the 
use  of  the  scalar  floor  a  in  the  definition  of  Q(f )  in  Eq.(4).  In  Fig.  16  and  17 
we  show  the  corresponding  results  for  a  60-kHz  bandwidth  pulse  centred  at  65 
kHz.  From  these  figures,  it  can  be  seen  that  the  agreement  of  the  output  pulses 
with  the  desired  pulses  is  very  good,  particularly  for  the  60-kHz  bandwidth. 

Instead  of  these  very  short  pulses,  we  can  apply  the  same  compensation  to 
produce  cosine-tapered  Chirp  pulses  of  1  millisecond  duration.  Figures  18  and 
19  show  the  resulting  Chirp  pulses  with  the  desired  Chirp  pulses  superimposed 
in  red  for  a  80-kHz  wide  Chirp  and  in  Figures  20  and  21  the  corresponding 
results  for  the  60-kHz  wide  Chirp.  Perhaps,  more  informative  than  these 
time  series  plots  are  the  spectrograms  obtained  by  using  short-time  Fourier 
Transforms  over  the  received  signal.  The  resulting  time/frequency  plots  are 
shown  in  Figs.  22  and  23.  As  can  be  seen  the  spectograms  of  the  recorded 
signals  exhibit  the  correct  frequency/time  variation.  The  second  arrival  which 
is  evident  in  the  plots  is  from  the  surface  reflection.  The  vertical  angle  for  this 
multipath  is  quite  large,  due  to  the  short-range  geometry  used  in  the  tank, 
and  since  the  frequency  response  of  the  transducer  at  this  vertical  angle  is 
quite  different  than  that  at  horizontal,  the  surface  reflection  is  not  simply  an 
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Figure  14:  The  time  series  resulting  from  the  input  waveform  of  Fig.  13  -  blue  is 
time  series,  red  is  desired  time  series.  The  desired  time  series  is  shown  at  twice 
the  hydrophone  sampling  rate. 


Figure  15:  The  measured  output  spectrum  corresponding  to  the 
pre-corn pensated  80-KHz  bandwidth  Sine  time  series. 
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Figure  16:  The  time  series  resulting  from  the  compensated  waveform  for  a 
60-kHz  bandwidth  Sine  pulse  -  blue  is  time  series,  red  is  desired  time  series.  The 
desired  time  series  is  shown  at  twice  the  hydrophone  sampling  rate. 


Figure  17:  The  spectrum  corresponding  to  the  output  60-KHz  bandwidth  Sine 
time  series. 
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Figure  18:  The  resulting  Chirp  time  series  (blue)  with  desired  series 
superimposed  in  red  -  80  kHz  bandwidth 


inverted  replica  of  the  direct  signal. 

It  is  also  possible  to  construct  a  pulse  sequence  by  amplitude  or  phase  mod¬ 
ulating  a  basic  pulse.  For  example,  in  Fig.  24  we  show  2  5-pulse  sequences 
(as  recorded  after  compensation)  constructed  from  a  Gaussian  pulse  (cr/  =  20 
kHz,  centre  frequency  =  65  kHz)  with  positive  and  negative  polarities.  The 
difference  between  the  2  input  pulse  sequences  used  for  Fig.  24  is  the  inter¬ 
pulse  spacing:  80/isecs  in  the  first  case  and  40/isecs  in  the  second  case.  For 
this  case  we  simply  created  the  compensated  time  series  by  superimposing 
the  compensated  waveforms  for  the  individual  pulses.  This  sequence  of  ±1 
follows  a  Barker  code  [3]  (1,1,1, -1,1).  In  Fig.  25  we  show  an  11-bit  Barker  Se¬ 
quence  (1,1,1,- 1,-1,- 1,1,- 1,-1, 1,-1)  constructed  from  a  sequence  of  4-cycle  62.5 
kHz  pulses.  The  spectrum  of  the  recorded  Barker  sequence  of  Fig.  25  is  shown 
in  Fig.  26.  In  Ref.  3  an  application  of  these  pulses  to  doppler  sonars  is  pre¬ 
sented. 


2.4  Some  Matched  Filtering  Results 

When  using  extended  pulses  such  as  the  one  millisecond  Chirps  or  the  mul¬ 
tipulse  Barker  Code  sequences,  it  is  anticipated  that  by  match-filtering  the 
time  series  the  original  extended  pulses  can  be  reduced  to  narrow  peaks  with 
respect  to  time  in  the  matched  filter  domain.  For  example,  in  Fig.  27  we  show 
the  matched  filter  results  (amplitude  of  the  complex  envelope  formed  using  the 
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Figure  19:  A  zoom  of  the  80-kHz  bandwidth  Chirp  time  series  (blue)  with 
desired  series  superimposed  in  red 


Figure  20:  The  resulting  Chirp  time  series  (blue)  with  desired  series 
superimposed  in  red  -  60  kHz  bandwidth 
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Figure  21:  A  zoom  of  the  60-kHz  bandwidth  Chirp  time  series  (blue)  with 
desired  series  superimposed  in  red 
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Figure  22:  Computed  time/frequency  plot  of  experimental  80-kHz  Chirp 
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Figure  23:  Computed  time/frequency  plot  of  experimental  60-kHz  Chirp 


matched  filter  time  series  and  its  Hilbert  transform)  when  we  cross-correlate 
the  recorded  time  series  with  the  theoretical  input  pulses  for  the  5  pulse  se¬ 
quence  of  Gaussian  pulses,  first  when  the  individual  pulses  are  separated  by 
80  ^seconds  and  by  40  ^seconds.  A  very  good  time  compression  has  resulted 
for  both  cases.  In  Fig.  28  we  show  the  results  of  match  filtering  the  recorded 
time  series  produced  with  the  11  bit  Barker  Code  62.5  kHz  cosine  pulse  with 
the  desired  pulse  sequence.  Once  again,  a  very  good  time  compression  has 
been  achieved.  In  these  matched  filtered  time  series  the  second  peak  corre¬ 
sponds  to  the  surface-reflected  arrival.  As  discussed  previously,  due  to  the 
vertical  azimuthal  dependence  of  the  transducer  frequency  response,  the  sur¬ 
face  reflection  does  not  cross-correlate  as  well  as  the  direct  arrival  with  the 
incident  pulse.  The  frequency  compensation  which  works  well  for  the  small 
vertical  angles  angles  is  not  correct  for  energy  at  higher  vertical  angles. 

As  another  example  of  this  we  consider  the  time  series  for  the  80  kHz  Chirp. 
In  this  case,  the  receiver  and  transducer  were  moved  to  approximately  3  feet 
below  the  water  surface.  In  this  case  the  pulses  from  the  direct  and  surface 
reflected  paths  overlap.  As  well,  there  are  other  reflections  from  the  back  and 
sides  of  the  tank.  Also,  a  small  aluminum  sphere  (air-filled)  (shown  in  Fig. 29) 
was  placed  behind  the  receiver.  The  echo  from  this  sphere  is  weak.  In  Fig. 30 
the  original  time  series  is  shown.  The  incident  and  surface  reflected  pulse 
overlap  with  each.  The  echo  from  the  sphere  can  not  be  identified. 

The  results  of  match-filtering  the  time  series  with  the  theoretical  Chirp  pulse 
is  shown  in  Fig.  31  (here  we  simply  show  the  matched  filter  results  and  not 
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Figure  24:  A  5  pulse  Barker  sequence  using  individual  Gaussian  pulses  with  a 
separation  of  20  samples  (80  psecs)  between  pulses  and  10  samples  (40  psecs) 
between  pulses) 
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Figure  26:  The  spectrum  of  the  Barker  code  time  series  of  Fig.  25 


20 


DRDC  Atlantic  TM  2005-022 


C  ross  correi  ati  on  C  ross  correl  ati  on 


1400 


Figure  27:  The  matched  filter  time  series  using  the  Barker  code  of  individual 
Gaussian  pulses  separated  by  80  psecs  and  40  psecs  respectively) 
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Figure  28:  The  matched  filter  time  series  from  the  11  Barker  code  using  a  62.5 
kHz  cosine  burst. 


Figure  29:  The  small  aluminum-shelled  sphere  that  was  deployed  in  tank 
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Figure  30:  The  timeseries  with  overlapping  incident  and  multipath  pulses. 


the  amplitude  of  the  complex  envelope)  and  a  zoom  of  the  region  near  the 
sphere  echo  is  shown.  Although,  the  echo  is  relatively  weak  it  is  certainly 
detectable  in  the  matched  filtered  time  series.  This  echo  sequence  contains 
the  specular  reflection  and  other  reflections,  possibly  from  the  flange  as  well 
as  circumferential  arrivals. 

The  matched  filtering  results  are  consistent  with  those  predicted  by  theory. 
For  example,  the  width  of  the  main  lobe  for  the  11-term  Barker  code  (using 
the  62.5  kHz  burst)  is  approximately  132  p  seconds  which  corresponds  to  an 
effective  bandwidth  of  approximately  15  kHz,  which  is  consistent  with  the 
spectrum  of  Fig.  26.  When  using  the  80-kHz  Chirp,  the  width  of  the  main 
matched-filter  lobe  is  40  p  seconds  corresponding  to  an  effective  bandwidth 
of  50  kHz,  which  is  also  consistent  with  the  spectrum  of  the  the  pulse  (even 
though  we  used  80  kHz  in  the  definition  of  the  Chirp,  the  effective  bandwidth 
is  less  due  to  the  properties  of  the  Chirp  and  the  cosine  weighting  applied). 
Also,  as  illustrated  with  the  echo  time  series  from  the  sphere,  the  use  of  the 
matched  filtering  has  improved  the  signal-to-noise  ratio  from  the  original  time 
series. 
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3  VERTICAL  BEAMPATTERN 


In  the  data  presented  so  far,  even  for  the  2  element  MMPP  arrays,  the  sonar 
was  deployed  with  the  axis  of  symmetry  vertically.  For  example,  for  a  single 
MMPP  the  middle  slot  would  be  oriented  in  a  horizontal  plane.  The  hor¬ 
izontal  beampattern  should  be  approximately  omnidirectional.  The  vertical 
beampattern,  however,  will  be  much  more  complicated.  To  measure  this  ver¬ 
tical  beampattern,  the  MMPP  was  held  horizontally  and  then  rotated  in  an 
uniform  fashion.  The  FFT  spectrum  of  the  recorded  signal  for  a  80-kHz  wide 
incident  pulse  (noncompensated)  is  shown  in  Fig. 33  as  a  function  of  the  ro¬ 
tation  angle.  It  is  clear  from  this  plot  there  is  significant  energy  from  the 
circular  endcaps  (approximately  60°  and  240°  in  the  plot)  of  the  transducer 
and  this  region  is  relatively  broad.  The  spectra  which  we  have  been  examining 
in  this  paper  correspond  to  the  broadside  direction  (approximately  150°  and 
330°).  It  can  be  seen  that  the  beamwidth  of  this  broadside  energy  is  quite 
narrow,  particularly  at  the  higher  frequencies.  This  would  mean  that  in  its 
vertical  orientation  there  would  be  an  omni-direction  horizontal  beam  and  a 
narrow  vertical  beamwidth.  On  the  other  hand,  if  the  transducer  is  oriented 
horizontally,  there  is  an  omnidirectional  vertical  beampattern  and  a  very  nar¬ 
row  horizontal  beam  pattern.  It  is  also  clear  that  the  endfire  directions  which 
would  have  a  fairly  broad  conical  beam  could  also  be  useful.  The  null  at  ap¬ 
proximately  30  kHz  for  the  broadside  direction  appears  to  be  somewhat  less 
for  the  endfire  direction. 

In  Fig.  34  we  use  the  compensated  80-kHz  Sine  waveform  and  consider  the 
vertical  beampattern  measurement  in  this  case.  The  compensation  has  been 
successful  for  the  broadside  direction  but  has  not  compensated  the  other  di¬ 
rections  correctly.  In  fact,  the  relative  amplitude  of  the  signal  from  the  endfire 
direction  appears  to  be  somewhat  higher  now  (as  compared  to  the  case  of  the 
uncompensated  spectra).  Finally,  in  Figs. 35  and  36  we  show  the  recorded  time 
series  and  corresponding  spectrum  for  a  compensated  Sine  pulse  (20-70  kHz) 
with  the  transducer’s  flat  circular  endcap  facing  the  hydrophone.  Here  the 
spectrum  S(f)  used  in  Eq.(5)  was  that  for  the  endfire  direction.  As  can  be 
seen,  the  compensation  for  the  endfire  orientation  has  been  very  successful. 
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Figure  33:  The  vertical  beampattern  of  the  transducer  for  an  uncompensated 
80-kHz  Sine  waveform 
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Figure  34:  The  vertical  beampattern  of  the  transducer  for  a  compensated 
80-kHz  Sine  waveform  for  the  broadside  direction 
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Figure  35:  Time  series  using  compensated  Sine  pulse  (20-70  kHz)  for  transducer 
endcap  facing  hydrophone 


Figure  36:  The  power  spectrum  of  the  resultant  signal 
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4  DISCUSSION  OF  RESULTS 


In  this  report  it  has  been  shown  the  high  frequency  MMPP  has  a  very  wideband 
response.  The  spectral  response  of  the  transducer  can  be  simply  compensated 
for  in  order  to  produce  specified  wideband  pulses,  either  very  impulsive  or 
extended  in  time  in  the  case  of  1  millisecond  Chirp  pulses.  Additionally,  pulse 
sequences  can  be  created  by  amplitude  modulating  individual  pulse  types.  As 
an  example  of  these  pulse  sequences,  we  used  Barker  code  sequences  for  the 
amplitudes.  It  was  shown  that  by  match  filtering  the  signals  recorded  at  the 
hydrophone  with  the  ideal  (before  compensation)  waveform,  either  a  pulse 
sequence  or  a  Chirp,  very  precise  arrival  times  (i.e. ,  very  time  compressed 
outputs)  with  small  sidelobe  levels  were  obtained.  Thus  the  high  frequency 
MMPP  shows  good  potential  for  the  production  of  specified  pulses  in  the 
frequency  band  35-115  kHz  (for  the  standard  vertical  orientation). 

The  vertical  beampattern  of  the  transducer  was  also  investigated.  It  was  found 
that  frequency/azimuthal  dependence  of  the  transducer  is  rather  complicated. 
In  this  report  we  have  emphasized  the  on-axis  response  of  the  transducer 
with  the  endcaps  oriented  vertically.  The  vertical  beampattern  measurements 
indicate  that  in  this  orientation  the  transducer  has  quite  a  narrow  vertical 
beampattern  and  an  almost  omni-directional  horizontal  beampattern  (from 
the  azimuthal  symmetry  of  the  transducer).  For  applications  where  a  broad 
vertical  beampattern  and  a  narrow  horizontal  beampattern  are  required,  these 
characteristics  can  be  obtained  by  orienting  the  transducer  with  the  endcaps 
horizontal  and  the  centre  slot  facing  the  hydrophone.  It  is  also  seen  from  the 
vertical  beampattern  measurements  that  the  compensating  spectrum  used  for 
the  on-axis  energy  is  not  correct  for  angles  sufficiently  away  from  this  direction, 
particularly  at  the  higher  frequencies.  There  are  advantages  to  orienting  the 
transducer  so  that  one  of  the  endcaps  is  facing  the  receiver.  This  orientation 
seems  to  provide  more  power  and  its  beampattern  is  conical  in  nature.  The  null 
at  30  kHz  is  relatively  less  than  for  the  broadside  orientation  and  a  20-70  kHz 
Sine  pulse  was  very  accurately  produced  for  this  orientation.  The  beampattern 
is  broader  for  this  orientation  and  thus  the  computed  compensating  spectrum 
should  be  valid  over  a  larger  range  of  angles. 

It  should  be  possible,  due  to  their  small  size,  to  construct  an  array  of  these 
high  frequency  MMPP  sources.  The  overall  frequency  and  beampattern  char¬ 
acteristics  of  the  resulting  system  will  depend  upon  the  coupling  between  the 
transducers,  the  wiring  (e.g.,  in  phase  or  out-of-phase)  of  the  transducers  and 
the  inter-element  spacing. 
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In  the  future  it  is  hoped  to  investigate  the  performance  of  all  sizes  of  the  MMPP 
in  longer-range  sea  experiments.  It  is  hoped  that  a  prototype  wideband  sonar 
system  consisting  of  possibly  an  array  of  sources  and  an  array  of  receiver 
elements  can  be  constructed. 
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